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ABSTRACT 


Sea surface temperature (SST) anomalies on the daily synoptic time 
scale from National Oceanic and Atmospheric Administration/Nationa| 
Environmental Satellite Service (NOAA/NESS) Global Operational SST Com- 
putation (GOSSTCOMP) data in the Eastern South Pacific are compared To 
deviations of similar time scale from the means of both the wind stress 
curl and the friction velocity cubed as derived from Fleet Numerical 
Weather Central's (FNWC) mercator global band surface wind field. 
Example patterns of those variables as well as mean and unaltered fields 
are presented here. High correlations between the sea surface tempera- 
ture deviations and the deviations of wind stress curl and friction 
velocity cubed were not attained at any log, owing largely to the 
paucity and poor quality of surtace wind and satellite data. Further 


investigation is needed. 
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|. INTRODUCTION 


The meteorological and oceanographic activity of The Southern Hemi- 
Sphere, although not nearty as widely studied as that of the Northern 
Hemisphere, is still of vital importance to the United States and the 
rest of the world. Its cause-and-effect impact on the weather and ocean 
situations of the North, as well as on vast resources there not well 
utilized, arouses scientific curiosity to understand the physical pro- 
cesses in the region. 

One facet of These activities which is of particular interest is sea 
Surface temperature (SST). This paper will attempt to determine if 
weekly averaged non-seasonal fluctuations of wind phenomena analyzed 
by Fleet Numerical Weather Central (FNWC) can be correlated with non- 
seasonal fluctuations of sea surface temperatures revealed by Global 
Operational SST Computation (GOSSTCOMP) data. The major eventual goal 
is To examine mechanisms responsible for synoptic scale fluctuations In 
SST. If changes in wind phenomena can be related to the changes in SST, 
Then temperature inferred from wind fields may serve as a substitute 
for satellite data and possibly for prediction purposes. 

The specific area studied for this report was the Eastern South 
meaeitic, Off the coasts of Ecuador, Peru, and Chile. 

The wind phenomena chosen for correlation were perturbations of the 
f——Enical component of the curl of the wind stress (curl (7) and friction 
velocity cubed ie Oe (Symbols used in this thesis are defined on the 
List of Symbols on page 9.) These parameters were used because the 


curl of Tau is responsible for Ekman convergence and divergence, and 


sage 





frictional velocity is related to mixing. Both of these can be found to 


Seni buipe FO sol Changes. 
Meridional Ekman advection (as represented by Tau.) could also be 


studied for relationship with SST, but it was not found To be practical 


with data used in this paper. 
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Pomeoe @eLOrRMENT OF BASTC THEORY 
INVOLVING U,%, CURL OF TAU, 


AND SST 


Starting with the first law of thermodynamics, 


: a + p24 =) where H is the rate of heating per unit mass. 


: at 
Assuming incompressibility, then 


da 


mo 


By expanding the total derivative of T and averaging it, while 
neglecting VeV5T and ee (horizontal derivatives), we get 
aT , oT , 38) ee L _ 

eet iS a5 (WIT!) = Q/p, Cy = RUZ), 

where R(z) is the absorption of solar radiation itn the sea, 


let H = oe and W'T! is proportional to the upward flux. 


lf we look at the mixed layer én ie, oe =O | = If. end sor0Z) 


dZ S 
becomes 
os. 2. (WITT) = R(z) with z> -h 
at OZ i 
Next integrate equation 3 itn depth from the surface to z = -h and 


neglect R(z), Then 


o (UO! 
= = [owtt) — «wim 
A, / ae dz LCwtT ds (WIT Joel which leads to 


fe) aT 22 _ 
oe f aa = ~(WIT')  +CWITT) 
O zi 


=) 


ot 


Then with strong winds it would be expected that TET (whieh 1S. pro= 


Beei!ional to U, multiplied by the air sea temperature difference) be 


larger and (WITT), Ciiiethsaenhonem | Ona! fo W) be more negative 


(Figure 1). Thus, enhanced upward flux at the surface and downward flux 


at -h would result in lowering the SST. 
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Figure |. Influence of wind stress (tT) on ocean temperature (71) 
Structure with depth. Solid line represents structure before U, 
changed mixed layer. Dashed line represents resu!ting temperature 
profile after being effected by strong wind stress. Note how the 
mixed layer has lowered the surface temperature and increased the 
depth of the mixed layer. 
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In the equation (2) if oa is not zero, there will be an effect of 


vertical advection. mE is the mean vertical Ekman velocity and is 
directly related to the vertical component of the curl ee 


ewe = —— curl. t 
3 pf Z 


In the Southern Hemisphere f is negative, and therefore, the current 
is directed 90° to the left of the wind. Also in the Southern Hemisphere 
where the winds in highs and lows rotate in opposite directions to that 
of the Northern Hemisphere, there is horizontal Ekman convergence and 
Sinking motion within a surface counter-clockwise circulation and hori- 
zontal Ekman divergence and rising motion within a surface clockwise cir- 
culation. With this convergence, associated warming should be expected 
when the sign of curl, t 1S positive. Conversely clockwise circulation 
should have divergence yielding cooling with curl of Tau being negative. 
This process of surface diverging (converging) producing cooler (warmer) 
water has been called Ekman pumping (Figure 2). 

The purpose of this thesis was to try to detect wind mixing (Ue) 


and/or Ekman pumping (curl, t) from FNWC wind data and the GOSSTCOMP data. 
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eer One 


Sea surface Temperature observations were obtained from National 
Oceanic and Atmospheric Administration/National Environmental Satellite 
Service (NOAA/NESS)'s GOSSTCOMP datum. 

The GOSSTCOMP data was a NOAA satellite scanning radiometer product 
which had been modified and composited in time. A detailed description 
of this product is found in NOAA/NESS Technical Memorandum #78, June 
'76. During the time frame of this report, past SST information (climato- 
logy) was used as one type of validity check. 

This data was supplied weekly by NOAA/NESS through Defense Mapping 
Agency--Inter-American Geodetic Survey (DMA-IAGS), Canal Zone, on com- 
Duter print-out paper. 

The area of interest was from the coast of South America to |00°West 
longitude, and from the Equator to 40° South latitude. SST was displayed 
in tenths of a degree every one-half degree of latitude and longitude 
elona with a confidence factor for that area. 

For the purpose of the study, SST's were picked manually every 
2 1/2° latitude and longitude in order to smooth the field, to relate 
it to wind data at every 2 |/2°, and, finally, to allow a reasonable 
time in card punching for computer graphics. Thus a SST field covering 
The area studied was produced for a period beginning 22 May 1974 and 
emaimg 26 May 1975, with a total of 535 maps. 

Marine surface wind observations came from FNWC global band mercator 


maps. This data covered the Northern Hemisphere up to 60° North and the 
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Southern Hemisphere to approximately 40° South. Every |2 hours wind 
fields of surface '"'u" component and '"'v" component of wind in centimeters 
per second for the time frame of the report were made available in tape 
form for this study courtesy of FNWC. (Thesis data available from 

J. B. Wickham, Department of Oceanography, NPS). 

The winds andtemperatures were defined at grid points which were 
equally spaced on the mercator projection and had to be converted to an 
equally spaced latitude-longitude (2 |/2°) grid. To do this, FNWC's 
Bessel interpolation scheme was used. 

For the purpose of comparison, mean monthly SST maps came from 
Fishing Information (1974-1975), a publication of National Marine 
Fisheries Service (NMFS). Furthermore, the NMFS also supplied data 


from the El Nino Watch Cruise (II February 1975 to 27 May 1975). 
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own Geoo I NoTOr TDN IA 


Using the wind field tapes supplied by FNWC, the wind stress (1) was 
calculated by the following equation, which involved two assumed constants, 
= —> a 
foe = tau = OC, lviv. 


For practical use this can be separated into components 


ee Cie os VANES ery 


= Z al 2 
ae Cy Gin Gey) ae 


where the density of air, p, , is .00122 gm/cm, and the drag coeffi- 


cient, Cy pls” ~OOl 5. 

The TY and by fields were calculated every I2 hours and Then averaged 
for seven days. The last date used in calculating the average corresponded 
TO a given SST data date. 

The main interest in this data was to be able to calculate the mean 


vertical component of the curl of Tau and mean te? It is belteved that 


the changes in SST can be related by equation 2, as noted in another 


secTion 
oT _ = oT 8 ee 
Ze a nee = aS Ce aD) 
C1) (2) (S) 
where term | is The change of mean temperature with time, term 2 is 


Meesor1iOne! to the curl of Tau, and term 3 is proportional to wind 


mixing. 


Ly 





The definition of friction velocity is 
> 
oa Uy, = dehy 1/2 (Glossary of Meteorology, 1959) 
a 


which leads To 


: A 2k 2 21/2 5/2 
9 i, (aa 
0 

a 


The finite difference form of the vertical component curl of Tau equation 


in spherical coordinates over a seven day average is 


nO. 
2 oe - | Copp get nae. VEIN) a SMO Te 
ij RIE), | Z2AX 2Ad 
R is the radius of earth in centimeters. 


AX and Ad equals the grid spacing which is 2 |/2° converted to radians. 
db is the latitude in radians. 


hy is the west to east component of wind stress averaged over a 
seven day period. 
iT is the south to north component of wind stress averaged over a 


seven day period. 


kK The vector k is a unit vector directed vertically. 


These equations were incorporated into an existing program for cal- 
culating mean tue and TIS. The results indicated that the curl Tau was 


of the magnitude lon? —— tor dyne/em, comparable to that of Hantel 
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(1970), which was 2.3 x |0 dyne/em”, and the cube of the friction 


velocity ee was of the magnitude of | x 10" (em/sec)?. 
The next step was to separate these seven day averaged fields into 


mean and perturbation values for the purpose of removing the seasonal 
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cycle. Using the basic equation A = A +A! where A is the mean and 
A' is the deviation of A from the mean, the weakly mean field was 
used to develop the annual mean field of each parameter (Figures 3, 4, 
5). Once mean fields were generated, they were subtracted from the 


original field as shown by equations below 


(with: i = number of days in a week and 
J = number of weeks in a year) 
, a | 7 
! = a ae 
(Curls eye Clea. az 2 S d (curl, 1) ,),] 
ceti@leet)'s= cunls tT = curl_ T 
id Z Ze 
Sf 7 
Br. ae | | 5 
a) =U, - sz ee Usa) J 
j=! i= | 
| Be 
Sy = sill las CSS Pca 
‘ie 
S51! =eSoe=-99 | 


The perturbation fields of the curl T, friction velocity cubed, 
and sea surface Temperature were put into contoured map form so that 
their spatial variations could be examined. From an examination of these 
maps a decision was made to concentrate the study on the mid-latitude 
regions. A line of data from 27.5° South to 39.5° South at 87.5° West 
was chosen because it was thought to be far enough from the coast To 
eliminate coastal effects and close enough at mid-latitude to be under 
the effect of mid-latitude storms. 

On this North-South line, consisting of five stations (Figure 6) 
Separated by 2 |/2° latitude, the seven-day average SST's were plotted 


for May 1974 to May 1975 (Figures 7A, 8A, 9A). A definite sinusoidal 


tek 
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Annual mean friction velocity cubed in 10° 
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shape was obvious in the data, especially the SST. Therefore it was 
necessary to remove the seasonal effect since the interest of this study 
was the synoptic time scale perturbations. 

Harmonic analysis of the data was done, and the mean and first 
three harmonics, eliminating the seasonal cycle, were removed, leaving 
the synoptic frequencies. This synoptic frequency part of the 
(curl, (Uy?) ', and SST' were then plotted in time series to verify 


that the seasonal effect had been removed (Figures 7B, 8B, 9B). 
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V. OBSERVATIONS 


A. WARM SST TONGUE 

A study of the sea surface temperature maps for the year period 
Showed a very apparent tongue of warm water migrating from the North to 
the South near the coasts of Ecuador and Peru. 

A moderately warm finger (17°C) began forming in September (Figures 
lO, I1, 12), sloping down from the northwest parallel to and about 5°S 
from the coast. It essentially remained stationary until November when 
it extended farther down the coast to about 20° South latitude. By the 
end of November the core was I8°C. At the end of December, the core 
warmed to 19°C and in January to 21°C. By February the maximum tempera- 
Ture was 23°C. The basic structure of this warm anomaly was still 
tongue-shaped from the North and extended to 20° South latitude. Its 
Temperature decreased in March, and as May developed, the tongue 
vanished totally. 

The change in temperature in this area (20°S, 75° West) was 8.4°C 
from 28 August 1974 (-3.8°C below mean) to 19 February 1975 (+4.6°C 
above mean). It would be difficult to attribute this SST difference 
wholly to seasonal change, or at least to the kind of seasonal change 
seen in neighboring areas which are notably smaller. 

B. COMPARISON OF NOAA/NESS GOSSTCOMP DATA WITH THAT 

OF NOAA/NMFS MEAN MONTHLY SST MAPS 

The products from NOAA/NESS GOSSTCOMP and NOAA/NMFS appeared some- 
what similar. The NMFS usual ly nd strong coastal gradients and North- 


7, 
South gradients at the equator. This was not shown on the GOSSTCOMP 
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Figure 10. Sea surface temperatures (°C). 
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Figure Il. Sea surface temperatures (°C). 
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data. Some significant anomalies were observed in May and November on 
both products. The basic trends in the GOSSTCOMP data were depicted in 
the NMFS generally, though not always consistently. 

It should be pointed out that the NMFS SST maps are means and did 
indicate that the data was sparse in this particular region. The 
GOSSTCOMP was an accumulated short time daily average. 

C. COMPARISONJOF SMOOTHED GOSSTCOMP DATA WITH DATA 

eo ele NINO WATCH CRUISE, || FEB 75 to 27 MAY 75 

This shipboard data was supplied by the NMFS from the University of 
Hawaii and the Scripps Institution of Oceanography investigation of 
El Nino (reference Wyrtki, et al. 1976). Data was extracted from this 
investigation that corresponded the closest In space and in time to 
available GOSSTCOMP data in The area of interest (Table |). 

Of the several random examples selected for comparison, the average 
absolute difference between the cruise data and GOSSTCOMP data was |.17°C. 
In NOAA/NESS TM 78, the indication for 0°-560° longitude from 20° South 
to 90° South was a mean absolute difference in ship measurement versus 
satellite measurements was 0.75°C. For the purpose of this study it was 
assumed that SST fluctuations equal to or greater than |°C in the 


GOSSTCOMP data were an accurate representation of the ocean conditions. 


Weeeoot’ AND (CURL OF TAU)! 
An observation of the maps of the deviation of the means indicated 
a relationship of (cur |, tT) to SST'. As previously mentioned in 
the Southern Hemisphere with associated positive(negative) values of the 
Ccurl | T)', Ekman convergence (divergence) leading to surface warm- 
ing (cooling) can be expected. Four cases of Ekman pumping were 


selected from the synoptic maps and summarized in Table ||. For example, 
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Figure 13. Sea surface temperature deviations from the mean in °C.j 
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Figure 14. Deviations from the mean of Viel Dee ond dyne/cm, 
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Figure 16. Deviations from the mean of curl, Tt in 107 dyne/em>. 
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Figure 17. Sea surface temperature deviations from the mean in °C. 


41 





95 W 92 W 85 W 62 W 73 W 78 W 





GUAYABUIL 


DEVIAI ICN OF CURL 


IZ55EB 75 





Figure 18. Deviations from the mean of curl, 16 ia fone dyne/em>. 
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Figure 19, Sea surface temperature deviations from the mean in °C. 
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Figure 20. Deviations from the mean of curl, eG lon dyne/em. 
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in Case |, at 30°S, 90°W from || September 1974, to 25 September 1974, 
the curl | T anomaly was strongly positive and the SST' increased during 
this time period from -5.1°C to -3.6°C. The other three cases show a 


Similar relationship. 


E. SST! AND (U,?)' 


5y4 


An examination of the anomaly maps indicated a relationship of (U, 
to SST'. As this wind increased above normal (or mean), lower values of 
SST' were seen. Two cases of mixing were selected from the synoptic maps 
and summarized in Table II]. For example, in Case | at 32.5°S, 75°W 
feces June TO 19 June 1974, the the anomaly was strong and the SST! 


decreased during this time period from +1.6°C to -2.7°C. The other case 


shows a similar relationship. 
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Figure 21. Sea surface temperature deviations from the mean in °C. 
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Figure 22, Deviations from the mean of ie in 107 cm/sec. 
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To establish a statistically significant relationship between synop- 
ile trequency fluctuations of SST', ent and (curl, T)' | some 
types of correlation was looked for at the five points of interest 
(Figure 6). The investigation points are at mid-latitudes and assumed 
to be unaffected by the coast. Al! are at 87.5° West longitude, at 
least 10° of longitude away from the coast. 

Available for use in this study were Biomedical Computer Programs 
from the Health Sciences Computing Facility, University of California, 
Los Angeles, California, at the Computer Facility, Naval Postgraduate 
Schoo!. From the P-series book, the stepwise regression program BMDPZR 
was chosen. All standard statistical parameters were available along 
with regression analysis and scatter plots. 

Of major interest was the normalized correlation matrix (values = 
O + 1!) which showed how well the synoptic frequency parts of (U4?) 
and (curl | t)' were related to synoptic frequency parts of SST'. 
Correlations r(t) were performed on these three variables with various 
lags of SST' with time. One lag equalled seven days. The equations 
used to correlate these variables were as follows: 

Memeemcool' (tf + lag) * (curl of t)'(7)) =t = time 
Beer (oSl'(t + lag) « Gow Cia) lag from 0 to + 4 weeks 

A lag of zero meant that the two variables were correlated at the 
Same time. Lags greater than (less than) zero meant that values of SST! 
were compared to values of ee and (curl tT)’ in the past (future). 


Inally, for each lag, r(t) was averaged over the year of data giving 
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time average lag-correlations (Table IV). A sample of the time averaged 
correlations at Stations #3 and #4 may be seen in Figures 27 and 28. 
Each variable was also correlated with itself to see how rapidly the 
parameter would change with time (Table V). Graphic representation of 
This auto-correlation shows that the synoptic frequency oscillations in 
3 


SU \! and Ceur 1s ie OCCU Reaten a SnoOtier time scale than that of 


Son (Figure 29). 
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VII. CONCLUSIONS 


The correlation values were not what were expected. The highest 
ear ion (t.9295) for (curt tT)’ occurred at a -2 week jag of 
SST! at point #4, and the highest correlation (+.3251) for ay 
Occurred at a -3 week lag of SST' at point #3; but these results can 
be given little weight since they are not part of a coherent pattern 
in the correlation matrix. 

The basic theory was that the curl of Tau when positive would pro- 
duce surface Ekman convergence of ocean waters with downwelling caus- 
ing warming (weak effect) of the sea surface. When the curl of Tau 
was negative, the surface Ekman divergence should come with upwelling 
and cooling (strong effect) of the sea surface. Thus the time average 
correlations between SST! and (curl t)' was expected to be posi- 
tive at positive lags. High surface friction wind (U,) should cause 
Surface mixing resulting in cooling. Thus SST! and Go should be 
negatively correlated at positive lag. 

In terms of having the "correct" sign of the correlations, station 
#3 came closest to the theoretical expectations in both the fields. 

If the data were perfect and could be completely trusted, it would 
appear that both processes, wind mixing and Ekman pumping, were occur- 
ring with approximately equal importance at the station. Qualitatively, 
about half the variance in non-seasonal SST fluctuations at two week 
Periods would be properly related to these two processes. Unfortunately 
the results are not similar at the other locations (Table IV). [In this 


Study therefore, this basic theory was not supported by results. 
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The best explanation for the discrepancy between data and theory was 
That the Southern Hemisphere lacks a stable data base. Figure 7B, SST 
deviations with the seasonal cycle removed, showed SST! amplitudes rang- 
ing from 0 to + 1.5°C. In Table 1, Section V.C, where GOSSTCOMP data was 
compared to ship board data, the average absolute error between that of 
cruise data and GOSSTCOMP data was |.17°C. Since the error was the same 
order of magnitude as the deviations, it was expected that these errors 
strongly influenced the SST deviations. There was no way in this report 
to verify FNWC's marine global band fields over the Pacific Ocean in the 
Southern Hemisphere. There were very few ocean observations in this 
region so there is a possibility of non-negligible errors in the surface 
wind fields as well. Thus, during the period of this study, the Southern 
Hemisphere products are of questionable accuracy. 

The number of ships in the South Pacific Ocean reporting SST values, 
a prime source for checking satellite observations in the GOSSTCOMP 
routine, was definitely low compared to that of the North Pacific. It 
was also pertinent that there was a scarcity of weather observation 
Platforms that fed data into FNWC from this very large region. Yet, even 
at present, these were the best available sources for study of the 
Phenomena in the Eastern Pacific Ocean in the Southern Hemisphere. 

The author strongly recommends that further investigation into this 
area be encouraged. More research in the Eastern South Pacific would 
Provide a broader study base to facilitate a better understanding of the 


Changes in this region's ocean temperatures. 
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